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Abstract 

We report on a search for electro-magnetic and/or hadronic showers (cascades) 
induced by high energy neutrinos in the data collected with the AMANDA II de- 
tector during the year 2000. The observed event rates are consistent with the ex- 
pectations for atmospheric neutrinos and muons. We place upper limits on a diffuse 
flux of extraterrestrial electron, tau and muon neutrinos. A flux of neutrinos with 
a spectrum $ oc E~ 2 which consists of an equal mix of all flavors, is limited to 
E 2 $(E) = 8.6 x 10~ 7 GeV cm~ 2 s _1 sr _1 at a 90% confidence level for a neutrino 
energy range 50 TeV to 5 PeV. We present bounds for specific extraterrestrial neu- 
trino flux predictions. Several of these models are ruled out. 

Key words: Neutrino Telescopes, Neutrino astronomy, AMANDA 
PACS: 95.55.Vj,95.85.Ry,96.40.Tv 
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1 Introduction 



The existence of high-energy extraterrestrial neutrinos is suggested by the 
observation of high-energy cosmic rays and gamma rays. Observation of neu- 
trinos could shed light on the production and acceleration mechanisms of 
cosmic-rays, which for energies above the "knee" (10 15 eV) remain not under- 
stood. Cosmic rays are thought to be accelerated at the shock fronts of galactic 
objects like supernova remnants, micro- quasars, and in extragalactic sources 
such as the cores and jets of active galactic nuclei (AGN) [1]. High energy 
protons accelerated in these objects may collide with the gas and radiation 
surrounding the acceleration region, or with matter or radiation between the 
source and the Earth. Charged pions, produced in the interaction, decay into 
highly energetic muon neutrinos and muons which further decay into electron 
neutrinos. Fermi acceleration of charged particles in magnetic shocks natu- 
rally leads to power-law spectra, E~ a , where a is typically close to 2. Hence, 
the spectrum of astrophysical neutrinos is harder than the spectrum of atmo- 
spheric neutrinos (~ E~ 3 ' 7 ) potentially allowing to distinguish the origin of 
the flux (see for example [2]). 

For a generic astrophysical neutrino source, one expects a ratio of neutrino 
fluxes $j, e : & u : § Ut ~ 1 : 2 : 0. Due to neutrino vacuum oscillations 
this ratio changes to & Ue : : $„ r « 1 : 1 : 1 by the time the neutrinos 
reach the Earth. Recently a search with the AMANDA detector was reported 
[3], resulting in the most restrictive upper limit on the diffuse flux of muon 
neutrinos (in the energy range 6 to 1000 TeV). Clearly, a high sensitivity to 
neutrinos of all neutrino flavors is desirable. The present paper reports on a 
search for a diffuse flux of neutrinos of all flavors performed using neutrino- 
induced cascades in AMANDA. 



2 The AMANDA Detector 



AMANDA-II [4] is a Cherenkov detector consisting of 677 photomultiplier 
tubes (PMTs) arranged on 19 strings. It is frozen into the Antarctic polar 
ice cap at a depth ranging mainly from 1500 to 2000 m. AMANDA detects 
high-energy neutrinos by observation of the Cherenkov light from charged 
particles produced in neutrino interactions. The detector was triggered when 
the number of PMTs with signal (hits) reaches 24 within a time-window of 
2.5 pis. 

The standard signatures are neutrino-induced muons from charged current 
(CC) interactions. The long range of high energy muons, which leads to 
large detectable signal event rates and good angular resolution results in re- 
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strictive bounds on neutrino point-sources [5]. 

Other signatures are hadronic and/or electro-magnetic cascades generated by 
CC interaction of u e and u T . Additional cascade events from all neutrino fla- 
vors are obtained from neutral current interactions. Good energy resolution, 
combined with low background from atmospheric neutrinos makes the study of 
cascades a feasible method to search for extraterrestrial high energy neutrinos. 



3 Update on Cascade Search with AMANDA-BIO 

Before the completion of AMANDA-II, the detector was operated in a smaller 
configuration. The results for the search of neutrino induced cascades in 130.1 
effective days of the 10-string AMANDA-BIO detector during 1997 have been 
reported before [6]. The same analysis has been applied to 221.1 effective days 
of experimental data collected during 1999. The AMANDA detector in 1999 
had three more strings than in 1997, yet data from these strings were not used 
in this analysis, so that the detector configuration used in the 1999 neutrino 
induced cascade search is very similar to that of 1997. 

Signal simulation for the analysis of 1999 data was improved to the standards 
reported in this letter. No events were found in the 1999 experimental data 
after all selection criteria had been applied. We will present results supposing 
a background of zero events. 

Using the procedure explained in this letter we obtain an upper limit on the 
number of signal events of /igo%=2.75 at a 90% confidence level, from which 
we calculate the limit on the flux of all neutrino flavors. Assuming a flux 
$ oc E~ 2 consisting of an equal mix of all flavors, one obtains an upper 
limit $90% = 8.9 x 10~ 6 GeVcm _2 s _1 sr _1 . In the calculation of this limit we 
included a systematic uncertainty on the signal detection efficiency of ±32%. 
About 90% of the simulated signal events for this limit have energies between 
5 and 300 TeV, while 5% have lower and 5% have higher energy. Differences 
between this result and the one obtained with 1997 experimental data [6] are 
due to the larger live-time in 1999 and improved simulation. 



4 Data Selection and Analysis for AMANDA-II 

The data set of the first year of AMANDA-II operation comprises 1.2 x 10 9 
triggered events collected over 238 days between February and November, 
2000, with 197 days live-time after correcting for detector dead-time. 
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The background of atmospheric muons was simulated with the air-shower 
simulation program CORSIKA (v5.7) [7] using the average winter air density 
at the South Pole and the QGSJET hadronic interaction model [8]. The cosmic 
ray composition was taken from [9]. All muons were propagated through the 
ice using the muon propagation program MMC (vl.0.5) [10]. The simulation of 
the detector response includes the propagation of Cherenkov photons through 
the ice as well as the response of the PMTs and the surface electronics. 

Besides generating unbiased background events, the simulation chain was op- 
timized to the higher energy threshold of this analysis. By demanding that 
atmospheric muons passing through the detector radiate a secondary with an 
energy of more than 3 TeV, the simulation speed is increased significantly. A 
sample equivalent to 920 days of atmospheric muon data was generated with 
the optimized simulation chain. 

The simulation of v e , and v T events was done using the signal genera- 
tion program ANIS (vl.O) [11]. The simulation includes CC and neutral cur- 
rent (NC) interactions as well as W~ production in the V e e~ channel near 
6.3 PeV (Glashow resonance). All relevant neutrino propagation effects inside 
the Earth, such as neutrino absorption or v T regeneration are included in the 
simulation. 

The data were reconstructed with methods described in Ref. [6]. Using the 
time information of all hits, a likelihood fit results in a vertex resolution for 
cascade-like events of about 5 m in the transverse coordinates (x,y) and slightly 
better in the depth coordinate (z). The reconstructed vertex position com- 
bined with a model for the energy dependent hit-pattern of cascades allows 
the reconstruction of the energy of the cascade using a likelihood method. 
The obtained energy resolution in log 10 .E lies between 0.1 and 0.2. The per- 
formance of the reconstruction methods have been verified using in situ light 
sources. 

Eight cuts were used to reduce the background from atmospheric muons by a 
factor ~ 10 9 . The different cuts are explained below. The cumulative fraction 
of events that passed the filter steps are summarized in Table 1. 

Since the energy spectrum of the background is falling steeply one obtains large 
systematic uncertainties from threshold effects in this analysis. For example, 
an uncertainty of ±30% in the photon detection efficiency translates to up to 
a factor 2 ±:L uncertainty in rate. Such effects can explain the discrepancies of 
Table 1 in passing efficiencies between atmospheric muon background simu- 
lation and experimental data. However, as will be shown later, the threshold 
effects are smaller for harder signal-like spectra. 

At the lowest filter levels (cuts 1 and 2), variables based on a rough first-guess 
vertex position reconstruction are used to reduce the number of background 
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# cut variable exp. MC 









atm. fj, 


atm. v e 


E- 2 u t 


1 


N ear l y /N hit < 0.05 


0.058 


0.033 


0.94 


0.63 


2 


N dir > 8 


0.030 


0.016 


0.89 


0.57 


3 


-^vertex < 7.1 


0.0027 


0.0012 


0.39 


0.35 


4 


-^energy ^ S . E reco 


0.0018 


0.00077 


0.35 


0.26 


5 


-60 < z rcco < 200 


0.0010 


5.9-10" 4 


0.28 


0.18 


6 


Preco^S ■ -^reco 


8.6-10" 4 


5.1-10- 4 


0.26 


0.15 


7 


L s > 0.94 


9.7-10- 6 


4.8-10- 6 


0.040 


0.091 


8 


E TCCO > 50 TeV 


8-10- 10 


7-10- 10 


2.8-10~ 5 


0.029 



Table 1 

Cumulative fraction of triggered events passing the cuts of this analysis. Values 
are given for experimental data, atmospheric muon background Monte Carlo (MC) 
simulation, atmospheric u e simulation and a v e signal simulation with an energy 
spectrum <I> oc E~ 2 . The flavor v e was chosen to illustrate the filter efficiencies, 
since interactions of v e always lead to cascade-like events. 

events by about a factor of 30. It is useful to define the time residual of a hit as 
the time delay of the hit time relative to the time expected from unscattered 
photons. The number of hits with a negative time residual, N eai \ y , divided 
by the number of all hits, AThits, m an event should be small. This first cut 
criterion is effective since early hits are not consistent with the expectation 
from cascades, while they are expected from long muon tracks. Cut 2 enforces 
that the number of so called direct hits, N^r (photons having a time residual 
between and 200 ns), is large. 

Cut 3 is a requirement on the reduced likelihood parameter resulting from 
the standard vertex fit, L VCTtcx < 7.1 (see also [6]). Note, that the likelihood 
parameter is, in analogy to a reduced x 2 , defined such that smaller values in- 
dicate a better fit result, hence a more signal-like event. In a similar manner, 
the resulting likelihood value from the energy fit, L energy, is used as a selection 
criterion (cut 4). However, since the average value of L energy has an energy 
dependence, the cut value is a function of the reconstructed energy, E reco . Cut 
5 on the reconstructed z coordinate, z TCCO , was introduced to remove events 
which are reconstructed outside AMANDA and in regions where the simu- 
lation of the ice properties for photon propagation is insufficient. While the 
upper boundary coincides roughly with the detector boundary, the lower value 
is about 100 m above the geometrical border of the detector. Restricting z reco 
improves significantly the description of the remaining experimental data (for 
example the reconstructed energy spectrum) [12]. Only events reconstructed 
with a radial distance to the detector z-axis, p rcco < 100 m, are accepted 
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Ap xy [ m ] L s 

Fig. 1. Normalized distribution of the three input variables L ver t eX ) cos 9^ and Ap xy 
as well as the resulting likelihood variable L s . Shown are experimental data as well 
as atmospheric muon and signal MC simulations after cut 6. 

(cut 6), unless their reconstructed energies lie above 10 TeV. For each decade 
in energy above 10 TeV one allows the maximal radial distance to grow by 
75 m. This reflects the fact that the cascade radius 1 , increases as a function 
of energy, while the expected amount of background decreases. 

Three discriminating variables are used to form the final quality parameter 
L s : 

1. The value of the reduced likelihood parameter resulting from the vertex fit, 
-^vertex- Note that this variable has been used previously in cut 3. 

2. The difference in the radial distance of the vertex position reconstructed 
with two different hit samples, Ap xy . While the first reconstruction is the 
regular vertex reconstruction using all hits, the second reconstruction uses 
only those hits outside a 60 m sphere around the vertex position resulting 
from the first reconstruction. Since the close-by hits typically contribute 
most to the likelihood function, their omission allows to test the stability 



We define the cascade event radius as the direction averaged distance from the 
vertex at which the average number of registered photon-electrons is equal to 1. 
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of the reconstruction result. If the underlying event is a neutrino-induced 
cascade, the second reconstruction results in a vertex position close to that 
of the first reconstruction. In case of a misidentified muon event, removing 
hits located close to the vertex typically results in a significantly different 
reconstructed position. 
3. The cosine of the angle of incidence cos 9^ as reconstructed with a muon- 
track fit. The muon-track fit assumes for the underlying likelihood parametriza- 
tion that the hit pattern originates from a long range muon track. The fit 
allows to reconstruct correctly a large fraction of the atmospheric muons. 

The final quality parameter is defined as a likelihood ratio: 

L = UiPijxi) (1) 



where % runs over the three variables. p h {h = s for signal and h = b for back- 
ground) are probability density functions defined &sp h (xi) = ft( x i)/(fi( x i) + 
fi(xi)). f h (xi) correspond to the probability density functions of the individual 
variables Xi for background due to atmospheric muons and signal consisting 
of a flux of v e with a spectral slope §(E U ) oc E~ 2 . They are obtained from 
simulations. 

The distributions of the individual variables as well as of the likelihood ratio L s 
are shown in Fig. 1 for experimental data, atmospheric muon background and 
signal simulations. The experimental distributions of Ap xy and L vertex approx- 
imately agree with those from the simulation while the distribution of cos 9^ 
shows some larger deviations. The deviation reflects an simplified description 
of the photon propagation through the dust layers in the ice [12]. The exper- 
imental L s distribution is not perfectly described by the atmospheric muon 
simulation, which is mainly related to the mis-match in the cos 9^ distribu- 
tion. The related uncertainties in the cut efficiencies are included in the final 
results. 

At this stage of the event selection one is left with events due to atmospheric 
muons, which happen to radiate (mostly through bremsstrahlung) a large frac- 
tion of their energy into a single electro-magnetic cascade. The reconstructed 
energy corresponds to that of the most energetic secondary-particle cascade 
produced in the near vicinity of the detector. To optimize the sensitivity of 
the analysis to an astrophysical flux of neutrinos, a further cut on the recon- 
structed energy, E reco , was introduced. 

The sensitivity is defined as the average upper limit on the neutrino flux 
obtained from a large number of identical experiments in the absence of signal 
[13,14]. The sensitivity was calculated for a flux of u e with spectrum oc E~ 2 . 
A flux of v e was used for optimization, since z/ e -induced events always have 



8 



4 ii i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 

Likelihood Parameter l_ s 

Fig. 2. Optimization of final cuts. The sensitivity for the diffuse flux of v e is shown 
as a function of cuts on -Ereco and L s . The coefficient next to the contour lines 
correspond to the average upper limit in units of (E/GeV)~ 2 • GeV -1 s _1 sr _1 cm -2 . 

cascade-like signatures. The sensitivity is shown in Fig. 2 as a function of 
the -E r eco and L s cut. L s > 0.94 and E reco > 50 TeV were chosen in this 
two dimensional optimization procedure such that the average upper limit 
is lowest. With these cuts the expected sensitivity for an E~ 2 spectrum of 
electron neutrinos is 4.6 x 10~ 7 (.E/GeV) -2 • GeV -1 s _1 sr _1 cm -2 . 



CM 




-1 

1 ^ 1 I 1 1 1 1 I 1 1 1 I I I 1 I 1 I 1 1 1 1 I 1 1 1 1 I 

4 4.5 5 5.5 6 6.5 

log 10 (E reco /GeV) 

Fig. 3. Distributions of reconstructed energies after all but the final energy cut. 
Shown are experimental data, atmospheric muon simulation and a hypothetical 
flux of astrophysical neutrinos. The final energy cut is indicated by the line with 
the arrow. 

The energy spectra of experimental data as well as signal and background 
simulations after all but the final energy cut are shown in Fig. 3. Note that 
the energy spectrum begins at 5 TeV, since this is the lowest energy for which 
the optimized background simulation is applicable. The number of events due 
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to simulated atmospheric muons was normalized to that observed in the ex- 
periment. 

One experimental event passes all cuts, while 0.961q;43 events are expected 
from atmospheric muons and a small contribution from atmospheric neutrinos. 

The spectrum as obtained from simulation of atmospheric muons passing cut 7 
was normalized to the number of experimental events resulting in an expecta- 
tion of 0.90tg;43 events due to atmospheric muons. The three main sources to 
the error are given by limited statistics of simulated atmospheric muon events 
(the error of was determined using the Feldman-Cousins method [14]), 
uncertainties in the cut efficiency (±20% obtained from variation of the cuts) 
and limited knowledge of the ice properties (±12% obtained from variation of 
the ice properties in the simulation). The total error was obtained by adding 
the individual errors in quadrature. 

The predicted event number from atmospheric neutrinos simulated according 
the flux of Lipari [15] is 0.06tg;o4 , where the uncertainties are mainly due 
to uncertainties in ice properties (error of ±0.03 obtained from variation of 
the ice properties in simulation), and in detection efficiencies of Cherenkov 
photons (lo!o2 obtained from variation of the photon detection sensitivity in 
the simulation). The theoretical uncertainties in the flux of atmospheric neu- 
trinos was estimated to be about 25% [16] and is small when compared with 
the other uncertainties. Again, the total error was obtained by adding the 
individual errors in quadrature. 

The uncertainty in the detection efficiency of neutrino events from an astro- 
physical flux with a spectral index a < 2 are estimated to be not larger than 
25 %. Because of the flatter energy spectrum, the uncertainties related to the 
energy threshold (such as the photon detection efficiency) result in smaller 
uncertainties in rate when compared to the uncertainties found for atmo- 
spheric neutrino events. The main sources of error are again uncertainties in 
the simulation of the ice properties (±15%) and the detection efficiencies of 
the Cherenkov photons (±20%). 

The experimental event which passed all selection criteria is shown in Fig. 4. 

The sensitivity of the detector to neutrinos can be characterized by its effective 
volume, Veg, or area, A e g, remaining after all cuts are applied. V e R represents 
the volume, in which neutrino interactions are observed with full efficiency 
while A e g represents the area with which a neutrino flux can be observed with 
full efficiency. While the concept of V c r is more intuitive because it relates to 
the geometrical size of the detector, the concept of A c r is more convenient for 
calculations of neutrino rates (see Eq. 2 in Sect. 5). 

Figure 5 shows V c q as obtained from simulation for all three neutrino flavors 
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Fig. 4. The experimental event which has passed all selection criteria is displayed 
from the side (left) and from above (right). Points represent PMTs, and shaded 
circles represent hit PMTs (early hits have darker shading, late hits have lighter 
shading). Larger circles represent larger registered amplitudes. The light pattern 
has the sphericity and time profile expected from a neutrino induced cascade. The 
arrow indicates the length scale. 




4 6 4 6 4 6 



log 10 (E v /GeV) 

Fig. 5. Effective volume for i/ e , j/„ and v T interactions as a function of the neutrino 
energy. The effective volume is shown without including Earth propagation effects 
(full line) and with Earth propagation effects (dashed line). 

as a function of the neutrino energy. The effective volume has been averaged 
over all neutrino arrival directions. As can be seen, V e g rises for energies above 
the threshold energy of 50 TeV. Above PeV-energies V e g decreases for v e and 
an effect related to both reduced filter efficiencies and neutrino absorption 
effects. In the case of v Ty the volume saturates because of regeneration effects: 
v T — ► r — ► v T and because of the event v T event topology: there is an increase in 
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log 10 (E v /GeV) 

Fig. 6. Effective area for u e , and v T interactions as a function of the neutrino 
energy. The effective area is shown without including Earth propagation effects 
(full line) and with Earth propagation effects (dashed line). 

detection probability for CC v T interactions (with energies above ~ 10 7 GeV) 
because the cascade from the hadronic vertex and the cascade arising from 
the subsequent tau decay are separated far enough in space to be detected 
independently. 

Fig. 6 shows A c ff as obtained from simulation for all three neutrino flavors as 
a function of the neutrino energy. Note that A c s is small because of the small 
neutrino interaction probability, which is included in the calculation of A c s 
(but not in V e s). 

The detector sensitivity varies only weakly as a function of the neutrino inci- 
dence angles. However, because of neutrino propagation effects effective area 
and volume are suppressed for neutrinos coming from positive declinations. 

The effect of the resonant increase of the cross-section for v e at the Glashow 
resonance is not included in A e g shown in Fig. 6. For energies between 10 6 7 
and 10 6 - 9 GeV the average effective area including Earth propagation effects 
is AZ% = 8.4 m 2 . 



5 Results 

Since no excess events have been observed above the expected backgrounds, 
upper limits on the flux of astrophysical neutrinos are calculated. The un- 
certainties in both background expectation and signal efficiency, as discussed 
above, are included in the calculation of the upper limits. We assume a mean 
background of 0.96 with a Gaussian distributed relative error of 73%, and an 
error on the signal detection efficiency of 25%. For a 90% confidence level an 
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Model 






v T 


V e + + V T 


^90% 
-^modcl 


1(T 6 x E- 2 


2.08 


0.811 


1.28 


4.18 


0.86 


SDSS [19] 


4.20 


1.91 


2.77 


8.88 


0.40 


SS Quasar [20] 


8.21 


3.57 


5.30 


17.08 


0.21 


SP u [21] 


33.0 


13.0 


20.5 


66.6 


0.054 


SP 1 [21] 


6.41 


2.34 


3.98 


12.7 


0.28 


P pp + pj [22] 


5.27 


1.57 


2.86 


9.70 


0.37 


P p 1 [22] 


0.84 


0.40 


0.56 


1.80 


1.99 


MPR [23] 


0.38 


0.18 


0.25 


0.81 


4.41 



Table 2 

Event rates and model rejection factors (MRF) for models of astrophysical neu- 
trino sources. The assumed upper limit on the number of signal events with all 
uncertainties incorporated is /Ugo% = 3.61 

upper limit on the number of signal events, /igo% — 3.61, is obtained using the 
Cousins-Highland [17] prescription implemented by Conrad et al. [18], with the 
unified Feldman-Cousins ordering [14]. Without any uncertainties the upper 
limit on the number of signal events would be 3.4. 

The effective area can be used to calculate the expected event numbers for 
any assumed flux of neutrinos of flavor i, <&i(E u ): 

AUdd = 4xttxT J2 f dE » HK)A^(E V ), (2) 

i=v e ,v,j„v T 



with T being the live-time. If iV mo d e i is larger than /zgo%, the model is ruled 
out at 90% CL. Table 2 summarizes the predicted event numbers for different 
models of hypothetical neutrino sources. Thereby, the spectral forms of z/ M 
and v e are assumed to be the same (the validity of this approximation is 
discussed in [12]). Furthermore, full mixing of neutrino flavors is assumed, 
hence <&„ e : $„ M : <&„ T = 1 : 1 : 1 as well as a ratio vjv — \. 

Electron neutrinos contribute about 50% to the total event rate, tau neutrinos 
about 30% and muon neutrinos about 20%. For the sum of all neutrino flavors 
the various predicted fluxes are shown in Fig. 7. 

The models by Stecker et al. [19] labeled "SDSS" and its update [20] "SS Q", 
as well as the models by Szabo and Protheroe [21] "SP u" and "SP 1" rep- 
resent models for neutrino production in the central region of Active Galac- 
tic Nuclei. As can be seen from Table 2, these models are ruled out with 
Jt 90% PS 0.05 — 0.4. Further shown are models for neutrino production in 

1 *modcl 

AGN jets: a calculation by Protheroe [22], which includes neutrino produc- 
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log 10 (E v /GeV) 



Fig. 7. Flux predictions for models of astrophysical neutrinos sources. Models rep- 
resented by dashed lines are excluded by the results of this work. Models fluxes 
represented by dotted lines are consistent with the experimental data. The labels 
are explained in the text. The solid line corresponds to the upper limit on a flux 
$ oc E~ 2 . 

tion through jry and pp collisions (models "P pp + pj" and "P P7") as well as 
an evaluation of the maximum flux due to a superposition of possible extra- 
galactic sources by Mannheim, Protheroe and Rachen [23] (model "MPR"). 
The latter two models are currently not excluded. 

For a neutrino flux of all flavors with spectrum oc E~ 2 one obtains the limit: 
£ 2 $ 90% = 8.6 x 10~ 7 GeV cm" 2 s" 1 sr" 1 . 

For such a spectrum, about 90% of the events detected have neutrino energies 
between 50 TeV and 5 PeV,with the remainder equally divided between the 
ranges above and below. The limit is shown in Fig. 7 as a solid line ranging 
from 50 TeV to 5 PeV. 

To illustrate the energy dependent sensitivity of the present analysis we restrict 
the energy range for integration of Eq. (2) to one decade. By assuming a 
benchmark flux $e (E) = $0 x (E/E )~ 2 x 9(0.5 - | hg(E/E )\) where $ = 
l/(GeV cm 2 s sr) represents the unit flux and G the Heaviside step-function 
(restricting the energy range to one decade), one obtains the number of events 
for a given central energy E : N event (E ) . The limiting flux at the energy E 
is then given by § 90% (E ) = % x fi 90% /N cvent (E ). The superposition of the 
limiting fluxes as a function of the central energy is shown in Fig. 8. For a flux 
$ oc E~ 2 the analysis has its largest sensitivity around 300 TeV. 

The mentioned strong increase in effective area at the energy of the Glashow 
resonance allows setting of a limit on the differential flux of V e at 6.3 PeV. 
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Fig. 8. Illustration of the energy dependency. The curved solid line represents a su- 
perposition of the limiting fluxes of a series of power law models <3? oc E~ 2 restricted 
to one decade in energy. The limit for a flux <E> oc E~ 2 without energy restriction is 
shown for comparison. The range of the dashed line represents the range of energies 
in which 90% of the signal events are detected while the embedded solid line rep- 
resents the energy range in which 50% of the signal events are detected, with the 
remainder equally divided between the ranges above and below. 



Re-optimizing the final energy cut for events interacting through the Glashow 
resonance results in an optimal cut, E rcco > 0.3 PeV. No experimental event 
has been observed in that energy range, which results in an upper limit on the 
number of signal events of /igo% = 2.65 assuming ±25 % uncertainties in the 
signal expectation and negligible background expectation. The limit on the 
flux at 6.3 PeV is: 

$ Pe (E = 6.3 PeV) = 5 x KT 20 GeV" 1 s" 1 sr' 1 cm" 2 . 



The transformation of this limit to a limit on a total neutrino flux is not 
completely straightforward, since the fraction of V e produced in the source is 
unknown. In cases of neutrino production through pp collisions one expects 
a ratio V e jv e 1. Hence, one would expect that about 1/6 of all neutrinos 
are V e . However, the V e produced in the pj — » nn + interaction through decay 
of the neutron carries only a very small fraction of the energy, and hence 
for most neutrino spectra contributes negligible to the high energy flux of 
neutrinos. For this flux of v e results mainly from neutrino oscillations. 

For maximal neutrino mixing, V e would constitute 1/9 of the total neutrino 
flux. If the mixing is non-maximal, that fraction would be smaller. 



1111111111111111111111111111111111111 
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6 Conclusion 



We have presented experimental limits on diffuse extragalactic neutrino fluxes. 
We find no evidence for neutrino-induced cascades above the backgrounds 
expected from atmospheric neutrinos and muons. In the energy range from 
50 TeV to 5 PeV, the presented limits on the diffuse flux are currently the 
most restrictive. We have compared our results to several model predictions 
for extragalactic neutrino fluxes and several of these models can be excluded. 

Results from the first phase of AMANDA, the 10-string sub-detector AMANDA- 
BIO, have been reported in [6] and an update to the analysis was presented 
above. Compared to AMANDA-BIO, the analysis presented here has a nearly 
ten times larger sensitivity, mainly achieved through using the larger volume 
of AMANDA-II and by extending the search to neutrinos from all neutrino 
directions. 

The limits presented here are also more than a factor of two below the AMANDA- 
BIO limit obtained by searching for neutrino-induced muons [3] and roughly 
as sensitive as the extension of that search using AMANDA-II 2000 data [24]. 
(Assuming a neutrino flavor ratio of 1:1:1, the numerical limits on the flux 
of neutrinos of a specific flavor (e.q. z/ M ) reported in the literature are 1/3 of 
the limits on the total flux of neutrinos.) The limits obtained from a search 
for cascade-like events by the Baikal collaboration [25] are about 50% less 
restrictive than the limits presented here. 

With the present analysis one obtains a large sensitivity to astrophysical neu- 
trinos of all flavors and in particular to electron and tau neutrinos. Hence, 
given the large sensitivity to muon neutrinos of other search channels, AMANDA 
can be considered an efficient all-flavor neutrino detector. 
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